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PRINCÍPIOS DE MÁXIMO EM EQUAÇÕES ÀS DERIVADAS 
PARCIAIS DE SEGUNDA ORDEM DE TIPO ELÍPTICO E SUA 
APLICAÇÃO EM TRANSMISSÃO DO CALOR (*) 


RESUMO 


Apresentam-se por via «fisica» os fundamentos dos 
principios de máximo em equações às derivadas parciais 
de segunda ordem de tipo elíptico e o enunciado de alguns 
teoremas importantes. Mostra-se como os teoremas são 
susceptiveis de aplicação no estabelecimento de condições 
de unicidade e na solução aproximada de eguações 
daquele tipo. 


0. INTRODUÇÃO 


A formulação matemática dum fonómeno físico 
envolve sempre aproximações que ficam subja- 
centes mesmo quando a solução matemática exacta 
se obtém nas premissas do modelo. Na maioria 
dos casos, porém, a solução é inacessível ou 
demasiado complexa e laboriosa para ter interesse 
prático, bastando o conhecimento duma solução 
matemática aproximada ou até, em muitas circuns- 
tâncias, apenas a forma geral da solução e algu- 
mas das suas propriedades. Em todos os casos é 
manifesto o interesse que haveria em obter uma 
estimativa do erro cometido. É nesta perspectiva 
que a presente contribuição se integra, e nela se 
pretende chamar a atenção para a aplicação prá- 
tica de alguns teoremas no domínio da teoria das 
equações às derivadas parciais de segunda ordem, 
de tipo elíptico e parabólico, que desempenham 


por J.J. DELGADO DOMINGOS 
Prof, 1.8.1. 
N.E.EN. do LAX., 


SYNOPSIS 


Using a physical reasoning the basic results for the 
maximum principles im partial differential equations 
of second order of elliptic type are established. À 
summary of some fundamental theorems is presented 
and it is shown how they cam be used for the demons- 
tration of conditions of unicity, and also on the cons- 
truction of apbbroximate solutions. 


um papel importante na demonstração de pro- 
priedades gerais, como a da unicidade da solução 
para certo tipo de condições impostas, e que, não 
apenas pela natureza dos textos em que se apre- 
sentam como pelo seu caracter muitas vezes 
abstracto, parecem ter eclipsado a possibilidade 
da sua utilização em situações bem concretas. 

Alguns dos teoremos que se vão apresentar 
têm aspectos que, sob o ponto de vista físico, são 
tão intuitivos que algumas conclusões, para uns 
quantos leitores, roçam pela banalidade. E muitas 
vezes terá sido esse carácter evidente, sob o 
ponto de vista físico, que inspirou a demonstra- 
ção matemática rigorosa, que por sua vez gerou 
demonstrações e conclusões que não possuem já, 
de modo algum, esse carácter evidente. 

Não é nosso objectivo fazer as demonstrações: 
apenas se enunciarão os teoremas e se farão 
algumas aplicações imediatas. 


(*) Matéria integrada num curso de pós-graduação regida no Imperial College (Universidade de Londres) a 


convite do Departamento de Engenharia Mecânica. 
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1 — RESULTADOS FISICAMENTE INTUITIVOS 


1.0 — Considere-se como exemplo introdutório 
a equação de Laplace 


1) v'6=0 


num domínio D, limitado. Um dos teóremas 
referidos, conhecido como «princípio de máximo» 
diz-nos que: 


«se & atingir o seu máximo, M, num ponto 
de D, então $==M em todo o domínio». 


O resultado é imediato se interpretarmos ú, por 
exemplo, como o campo das temperaturas num 
sólido homogéneo e isótropo de condutibilidade 
térmica constante e sem fontes de calor. Efecti- 
vamente, em regime estacionário, o campo das 
temperaturas ( rege-se por 1). 

Ora se o regime é estacionário, o fluxo de calor 


= — k grad q 


não varia no tempo; se 6 tem um máximo em 
M, pelo segundo princípio da Termodinâmica 
haverá dos pontos a temperatura mais elevada 
um fluxo de calor para os restantes, a tempera- 
tura inferior. Como em D não há fontes de 
calor, pelo primeiro princípio da Termodinâmica 
a temperatura em M diminui, o que contradiz a 
hipótese base do regime ser estacionário, único 
caso em que 1) representa o campo das tempe- 
raturas. 

Conclui-se pois que q=k grad $==0 e por- 
tanto é só pode ser constante. 

O raciocínio repete-se «mutatis mutandis» para 
o caso dum mínimo. 

Se é no interior do domínio não é constante, 
como também não pode possuir em D nem 
máximo nem mínimo, conclui-se que estes se 
verificarão na sua fronteira B (o que é também 
fisicamente intuitivo). 


A 


1.2 O teorema é susceptivel duma aplicação 
prática imediata do seguinte modo: 


Suponhamos que se pretende determinar é 
tal que 


eu U 
x GB 


1) V$/=0 
6 =F(x) 
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E que se conhece uma função O, obtida por 
exemplo da solução dum problema análogo, 
em que 


2) vV 0-0 x CD 


O=F(x) xcCB 

Se F, + E,, O não é obviamente solução de 1). 
Todavia, o teorema permite-nos conhecer o erro 
máximo que resulta de adoptar (O) como solução 
de 1). 


Efectivamente, seja 


v=6—6) 
pelo que vv = () 


Como o máximo e o mímino de v se verificam 
na fronteira B, será 


Vmaz = max (F — F) 
Vmin==min(F —F) 


e como v representa o erro de tomar O como 
solução de 1) as expressões anteriores dão os 
limites do erro máximo 


min(F, —F) < cuax SS max(F — EF) (*) 


max 


1.2.1 —- O teorema é ainda usado implicita- 
mente sempre que se exprime matematicamente 
a forma do domínio ou a condição imposta: 
quando se considera p. e. uma placa como limi- 
tada por dois planos a situação física é mais 
complexa porque as superfícies nunca são rigo- 
rosamente planas, no sentido matemático; de 
igual modo, a temperatura, mesmo admitida 
uniforme, nunca o é rigorosamente, de que resul- 
tará que em relação à solução matemática obtida 
existirá sempre um certo erro. O teorema permi- 
te-nos fixar os limites do erro, do seguinte modo; 


(1) Obviamente, se min ( F;, —-F)< 0, o maior desvio 
por defeito será min (F|/—F) e por excesso será max 
(F|— F) pois que na expressão se consideram os sinais 
respectivos. Considerando | F|—F | será 


imax = E max(F, — E) 
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sejam A e A” as superfícies reais que definem um 
domínio D' de Fronteira B, +B, (Fig. 1). 


Fig. 1 


Considerem-se agora os planos C e C” defi- 
nindo um domínio D, de fronteira B. e tal 
que D, € D”, e análogamento os planos E e E” 
definindo Dp, de fronteira B.etalque D'CD,.- 

Admitamos que a condição imposta em B, 
é T=T, e em B,,T=T, de que resultaria a 
solução T==F(P). 

Admitamos agora as mesmas condiçõesimpostas, 
respectivamente em B, e B.,eseja T, —F (P) 
a solução, e analogamente em By, e By, € seja 
T, = F, (P) a solução nesse caso 

Do teorema resulta que 


pelo que o erro máximo de ter substituído a 
superfície real por planos seria 
Cams E max(T —T,) 

De modo inteiramente análogo se consideraria 
o caso se ter tomado como condição imposta 
uma temperatura uniforme quando a não unifor- 
midade é E ==+ |AT|. A conclusão seria que 
o erro máximo na solução seria menor ou 
igual a £,. 

O erro máximo resultante da presença simul- 
tânea de irregularidade da superfície e da não 
uniformidade da temperatura obtém-se imediata- 
mente por combinação dos anteriores. 
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1.2.2 —- O resultado anterior mostra como se 
poderá formalmente demonstrar que a solução 
depende continuamente dos dados do problema, 
propriedade fundamental na descrição matemática 
de qualquer sistema físico. 


1.2.3— Como já se referiu, o princípio do 
máximo é utilizado muitas vezes na demonstração 
na unicidade da solução. A título exemplificativo 
considere-se o caso mais simples de condições 
de DIRICHELET — Temperatura imposta em toda 
a fronteira B do domínio (*). 

O problema formula-se matemáticamente como: 

Determinar & tal que 


3) vV6=0 xCD 


4) 6=F(x) x CB 

Admitamos que %, e &, são duas funções satis- 
fazendo 3) e 4). 

Então, dada a linearidade de 3) 


51 w= 0, — PD 


também satisfaz yº v=0 
Por outro lado, de 5) e 4) 


v=0 


Pelo teorema enunciado o máximo e o minimo 
de v verificam-se na fronteira; como v=-0 na 
fronteira, segue-se que v é idênticamente nulo 
e portanto 


em todo o domínio o que prova a unicidade da 
solução para condições de Dirichelet. 

Para condições de segunda e terceira espécie 
a demonstração é análoga a partir da generali- 
zação do teorema, como adiante se referira. 


() A demonstração matemática rigorosa implica 
algumas restrições quanto à regularidade do domínio e 
quanto à continuidade e diferenciabilidade da solução. 
Fisicamente, pelo menos no dominio da transmissão do 
calor, tais condições verificam-se sempre e não se analisam, 
dada a finalidade da exposição em que não cabe a conside- 


ração de soluções generalizadas. 
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1.3 Até agora limitamo-nos a considerar a 
equação de Laplace. Sucede que, quando se 
procuram soluções aproximadas para 


6) v$=0 


sujeito a certos tipos de condições aos limites, 
se é levado a propor como solução aproximada 


P- 4 
mas em que: 


7) ve=R(G)£o 


pelo que 2 não verifica a equação diferencial, 
dada a existência do resíduo 


R==R (x) 


Por hipótese suponhamos que v verifica as 
condições impostas na fronteira e que R(x)< O 
em todo o domínio. 

Atendendo ao significado físico de 7), 2 corres- 
ponde ao campo das temperaturas num sólido 
com fontes de calor de potência 


8) Q=-—çeR 


se R< 0 Q>0, pelo que fisicamente se 
conclui que 
> 6 


pois que é corresponde ao campo das tempera- 
turas num sólido idêntico, sujeito aos mesmos 
constrangimentos na fronteira, mas sem fontes 
de calor positivas. Estas tornam a temperatura 
necessáriamente superior em cada ponto. 

Análogamente, se R > 0, as fontes serão nega- 
tivos (poços) e 

0<y 


Com o exemplo anterior deu-se sentido físico 
ao resíduo R e introduziu-se a equação de Poisson 


9 
A dined 
Naturalmente, se as condições fronteiras forem 


idênticas entre 7 que satisfaz a equação anterior 
e | que satisfaz 


vy=Ri 
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existirão as seguintes relações: 


Se R,< R 
R,>2R 


>? 
LD 


o que é fisicamente evidente em face do signifi- 
cado de R, p. ex. em condução do calor. 


1.3.1 — Interpretando R como uma fonte de 
calor, os seguintes princípios de máximo são 
também evidentes: 


Se R>0, o máximo de 7 verifica-se na fron- 
teira (pois as fontes são negativas) 


Se R< 0, o mínimo de ? verifica-se na fron- 
teira. 


2.— GENERALIZAÇÃO DO PRINCÍPIO DO 
MÁXIMO PARA EQUAÇÕES DE 
SEGUNDA ORDEM DO TIPO ELÍPTICO 


2.0 — Introduzido heuristicamente o princípio 
do máximo para a equação de Laplace e esboçado 
como se constroem limites para a sua solução 
e para a equação de Poisson, apresentar-se-á 
seguidamente a generalização dos resultados para 
equações mais gerais de segunda ordem do 
tipo elíptico (º). 


2.1 — Considerem-se os operadores diferenciais 
de segunda ordem do tipo 


à 


n 
BÃO seo DES ms 
» O n) dx, 0x; 


ij=1 
e façamos 
1 
e id 


pelo que a expressão anterior se pode escrever 


n d 


o" 
L= 3 BR eras) ——— 


ij=1 dado 


(2.1) 


ã., = à, 


(') Não se apresentam as demonstrações porque a fina- 
lidade é apresentar teoremas susceptiveis de aplicação 
prática e não o desenvolvimento do formalismo que per- 
mitiu a sua dedução. 
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o operador (2.1) diz-se elíptico no ponto 
Ga O SR Xq) se a matriz simétrica 
dg * CM to A 
mo. 
Aly) = | * hj 
Sa A 8 dr Mia 


for positiva definida no ponto x. 
O operador é elíptico no domínio D se for 
elíptico em cada ponto do domínio. Tal significa que 


sendo é,, É, qualquer par de números reais. 


O operador é uniformemente elíptico em D 
se existir uma constante positiva |, tal que 


» (x) >, para todoo xCD 


2.2 — Considere-se a desigualdade 


od“u 
Za L jul = a. 
(2.2) [ul - ad 
du; 
ps 2 b, dx so (8) 2 A 
| 
se u for estacionário num ponto z,...z, será 


nesse ponto 


du, 
(5) = () j==1...fi 


e se esse ponto for um máximo deverá ser 


Desigualdade que não pode ser satisfeita no 
domínio se o operador for elíptico nesse domínio. 
Deste modo, qualquer função que verifique (2.2) 
não pode ter um máximo no dominio. 
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Admitiremos em toda a exposição seguinte 
que u é uma função contínua possuindo derivadas 
contínuas até pelo menos à segunda ordem 
(u é de classe C(?), As funções que figuram no 
operador são contínuas, o domínio D é limitado 
e a sua fronteira B é contínua, pelo menos 


por troços. (*) 
Teorema 1: seja 


du 


ja > 607020 


2 
Eds Taty" 


0x; dx; 


num domínio D em que L é uniformemente elíp- 
tico. Supunhamos que os coeficientes a, e b; são 
uniformemente limitados. Então, se u atingir um 
máximo M num ponto de D,u=M em D [2] 


Teorema 2: seja 


q? 
(LH hQg) mi = > ai; ru 


il dx; 


pe 2 b, (x) ei + h(x)u > 0 


com h<. 0, L uniformemente elíptico em D e 
aj, b; h limitados. Se u tiver um máximo não 
negativo no interior de D, então u=M [2] 


Teorema 3: seja 


(L4-h) [u]=£ xCD 
UU palx)u=8, x CB, 
dm 
U=2, x Cb, 
B=B,+B, 


Se L é uniformemente elíptico, h (x) <. O é limi- 
tado, e «(x) > 0 o problema tem uma solução 
única excepto quando h=z=0 e B,=0, caso 
em que a solução é única a menos numa cons- 
tante [3] 


(1) No caso do teorema 3 exige-se que em cada ponto 
da fronteira exista uma hiperesfera tangente totalmente 
contida em D. 

Sob o aspecto aplicado, as restrições não são impor- 
tantes. Sob o aspecto formal, algumas restrições podem 
remover-se. 
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Teorema 4: seja u tal que: 


(L+-h) [u]=f (x) xCD 
E patju— 8, (x) xCB, 
u= go (x) xCB, e B=B,+B, 


e z, tal que: 


(L+A) [z,] <£ (x) xCD 
dz, 
sn + <(x)z > gy (x) xCB, 
Z, > Bo(x) x CB, 
B=B, + B, 


então, seh (x) < Oea(x) >0 (*) 


u(x)< az (x) xCD 
Analogamente, se z, é tal que 
(L + h) [2,5] > fix xCD 


0 Z, 
E id «a l)z<g xCB, 


Z; S 8 x Cb, 


eh<0,4«<o 
então Za (x) <u (x) 


Teorema 5: seja À um valor próprio de 


du 
C++ [= Da, 6) ——— + 
Vi dx, dx; 
du 
+2b6 + Iht)+ikG|u=o 
i dx, xCD 
du 
——— “= (x) = 0) XcB, 
dn 
==) x CB, 


(') A condição Z — O verifica-se sempre em trans- 
missão do calor como consequência do segundo princípio 
da Termodinâmica. 
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B+B,=B 


ek (x)>0,n normal exterior. 
Seja w (x) numa função positiva em D + Be 
tal que 


da +alx)wz0 xCB, 
9 mn 
Então, a parte real de À 
Ro)> inf —L+th)w 
xcD kw 


(isto é, a parte real do valor próprio é superior 


ho 
= rh calculado para 


ao menor valor de - 
W 


todos os pontos do domínio). 


3. LIMITES PARA AS DERIVADAS DE EFUN- 
CÕES HARMÔNICAS 


3.0. Qualquer função u que satisfaça a equação 
de Laplace <7'u =0 no domínio D é, por 
definição, harmônica nesse dominio. 

De igual modo, se a função u verifica o" uz0 
em D é subharmónica em D, ou simplesmente 
subharmónica. Se «7º u < O será, pelo contrário 
superharmónica, 

Do teorema do máximo resulta imediatamente 
que se v é harmónica, v“ é subharmónica. Efec- 
tivamente 


Tv (v)=2vvyv+ 2 24 (50) 
i i 


4 / dv 
=13(5) co 
i dxi 


o que significa, portanto, que 


WE % “DD 
se v=-0 xe Bs 
ouv= il | CB 


3.1—Se u é harmónica, também o é cada uma 


; du 

das suas derivadas v, = «— 
A OX; 

se verifica, pelo que às suas derivadas também 
é aplicável o princípio do máximo. Em particular, 


, como facilmente 
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como à soma de funções subharmónicas ainda é 
subharmónica, segue-se que 


(8) 


é também subharmónica. 
Em condução do calor, como a densidade de 
fluxo é 


E rad u 


q=-—kgrad T 


e em regime estacionário sem fontes T é harmó- 
nico, q é superharmónico, pois k => 0. 


3.2.— Seja u (x,y) uma função harmónica em 
D, e Kp um circulo de centro (x,, y,) e raio r 
contido em D. 


Então 
Jlvuaxay= [| aiy (grad u) dx dy = 
K, K, 
A q 
CR OT 


sendo ds o elemento de arco ao longo de c,, 
circunferência que limita K, . 
Em coordenadas polares 


7) “du 
pas -[. a 


q 


como r O 


Integrando agora entre o e R e permutando a 
ordem de integração 


aa Es 
== / uR,)Rdi-zru(x,, y)=0 


o 
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ou seja 
1 dev 
u(xo, Yo) = u(R,0) Rdó = 
Ko 1 Yo TA 
o 
seje + uds 
27R c, 


expressão que constitui o teorema do valor 
médio para funções harmónicas. 

O resultado pode ainda exprimir-se do seguinte 
modo, multiplicando por R integrando de novo 
entre O e R: 


udA 


u(x,, Yo= 


) | 
considerando agora que pi também é uma fun- 
x 


ção harmónica resulta que 
5) e / 
=—— = du 
ç o am 2 gm] FT ada 


e pelo teorema da divergência aplicado ao 
segundo membro 


du 1 ar | 
dx)X Yo <R u (x, +R cos$,y, + 
o +Rsenf) Rcos9 d9 


Como se verifica o princípio de máximo, pode 
afirmar-se que 


m <=. u<.M 
sendo m e M respectivamente o minimo e o 
máximo de u em B. Por consequência a desi- 


gualdade verifica-se em D e portanto em K,. Por 
outro lado, se u é harmônica, 


j 
u + OM 
também o é pelo que: 


Den ca f e-tore) 
dy) XuVTo FKIoL 2 


Reosfdó< + 
RR 


o 


cos9d6 


u— É IM dem) 
2 


171 


e como u< M 


tai A Mm cui cos 9 d 6 = 
CF Gr TR 2 o 
2 (M-m) 
TR 


Dado que a direcção x escolhida é arbitrária, 
o limite ainda é válido para o caso em que a 
direcção escolhida coincide com o gradiente de u, 
pelo que 
(M-m) 


| 
| grad u (x, o — 
8 y) E 


sendo R a distância mínima do ponto (x,y) à 
fronteira de D, 


3.3 — No caso dum domínio tridimensional 
obter-se-ia 


> a 
radu (x,v,z) |< — (M-—m 
E (x, y42) | < (Mm) 


Sendo R o raio da maior esfera de centro em 
(x, y,Z) contida no domínio e Me m o máximo 
e o mínimo de u em B [3). 


3.3.1 — O limite indicado em 3.2 e 3.3 perde 
o seu valor à medida que o ponto se aproxima 
da fronteira porque então R — 0. 

Suponhamos, todavia, que a condição imposta 
na fronteira é 


u=g (x,y) 


e que g (x,y) é de classe C“em D+B. 

Admitimos ainda que a fronteira não possui 
ângulos reinterantes de tal modo que por cada 
x C B é possível fazer passar um círculo de raio 
po exterior a D. Considere-se ainda, concêntrico 
com o círculo anterior, um círculo de raio R tal 
que nele se contenha todo o domínio D. 

Demonstra-se então /3] que para qualquer 
ponto : 


erad u(x, y)|S cmo | grade | + + pos p [78] E. 


em 
Rê a 


O Cn 


4 R 
aplog ( ) 
Ê 


og R 


ai 
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Em particular, se o domínio for convexo, ? 
pode tender para coe obtém-se 


grad u | = max 
em B 


d 
ad + — sup q 
sra | pa À cá 


válido para qualquer número de dimensões, e 
sendo d a distância máxima entre dois pontos 
quaisquer em D., 


4 — ALGUMAS APLICAÇÕES 


4.0— O teorema 3 tem interesse imediato na 
medida em que nos garante a unicidade da solu- 
ção para as condições fronteiras habituais em 
transmissão do calor. 

Os teoremas 1, 2,e 4 permitem não só justi- 
ficar, como determinar os limites dos erros que 
se cometem nas simplificações quanto a geo- 
metria e condições impostas que é usual fazer 
nas aplicações. Em particular, a possibilidade de 
estimativa do erro permitirá a introdução segura 
de outras aproximações que a natureza e finali- 
dade da solução aconselham sem a incerteza que 
por vezes decorre da mera intuição. 

Uma aplicação imediata do teorema 5 encontra- 
-se em problemas de regime transiente em con- 
dução do calor em que a solução se procura pelo 
método de separação de variáveis. O teorema 
permite uma estimativa do limite inferior da 
velocidade de arrefecimento. Combinado com o 
método de Ritz, que forma um limite superior, 
possibilita o estabelecimento expedito de limites 
superior e inferior ; limites do mesmo tipo podem 
ser determinados para a extensão da região de 
entrada térmica numa conduta de secção arbi- 
trária. Citam-se dois exemplos em transmissão 
do calor, que se poderiam multiplicar e estender 
a outros domínios. 

41 — A título de exemplo, considere-se a deter- 
minação expedita dos limites de temperatura 
numa placa rectangular em que duas faces estão 
à temperatura T=-0 e as outrasa T==1. (Fig. 2) 
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UMA DAS DUAS PONTES ROLANTES DE 140t - 30t/20,5 FORNECIDAS À COMISION 
FEDERAL DE EITCTRICIDAD, S.A. PARA A SUA CENTRAL DE MALPASO (MÉXICO) 


PONTES ROLANTES, GUINDASTES E 
APAR. DE ELEVAÇÃO ESPECIAIS Projecto e fabrico 


TURBINAS HIDRÁULICAS ——— Fabrico segundo licença de A.C. M. de Vevey, S.A. 
TURBINAS A VAPOR ——— — Fabrico segundo licença de Brown Boveri, Cie. 
CALDEIRAS A VAPOR —————— Projecto e fabrico segundo licença de Foster 
| | Wheeler, Co. 
EQUIPAMENTOS E INSTALAÇÕES 
INDUSTRIAIS 


CONSTRUÇÕES METALOMECANICAS MAGQUE :... 


ALVERCA DO RIBATEJO- PORTUGAL 
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à ABES SR SR Sa 


Y 


SACOR — REFINARIA DO NORTE 
TANQUES PARA LUBRIFICANTES 


ALTERNADORES (de potência superior a 2000 kVA), COMPENSADORES 
SÍNCRONOS, MOTORES DE GRANDE POTÉNCIA . (fabricação parcial). 
EQUIPAMENTOS PARA INDÚSTRIAS QUÍMICAS E PETROQUÍMICAS 


Tubagens, tremonhas, reactores, permutadores, cisternas, etc. 


EQUIPAMENTOS PARA INDÚSTRIAS PETROLÍFERAS 


Torres, balões, permutadores, filtros, aerorefrigerantes, tubagens, etc. 


EQUIPAMENTOS PARA INDUSTRIAS CIMENTEIRAS 

Fornos, moinhos, separadores, arrefecedores, etc. 

RESERVATÓRIOS 

Le tecto fixo e tecto flutnante, de alta e baixa pressão. Para gases liquefeitos. Gasúmetros 


EQUIPAMENTOS DE MANUTENÇÃO E TRANSPORTE 


blevadores de alcatruzes, elevadores de palhetas, transportadores de tela metílica e de 


borracha, contentores, etc. 
ESTRUTURAS METALICAS 
Pontes, cimbres, coberturas, pisos, edifícios industriais, etc. 


MONTAGENS PESADAS 


SEPS SOCIEDADE DE CONSTRUÇÕES ELECTRO-MECÂNICAS, S.A.R.L. 
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APARTADO 8 * LEÇA do BAILIO «TELEX SEPSA P 2616 


O campo das temperaturas satistará 
4.1) vºT=0 
e às condições aos limites indicadas 


4.1.1 — Em primeira aproximação, vejamos o 
efeito da placa ser finita segundo uma das 
dimensões. 


Se a —º2, o problema seria unidimensional] 
e teríamos 


T, )=1— 
18 b 
Vejamos a relação entre Te T,: 
como JºT,=o0 
e T, satisfaz as condições fronteiras em y =0 e 
y==b, o desvio máximo verificar-se-á em x==0 e 
X= a 


Como imediatamente se conclui, 


T<Temx=0 
E | 


x== a 
pelo que T, não é limite superior nem inferior 
de T. 

Anãlogamente se conclui para T, se T, for 
solução do problema unidimensional 


Rede à 

dx? 

[=1 yv=o0 
T=0 y==a 


4.1.2 — Consideremos agora a solução 
Tt=T.T, 


Como T, = T,(y)e T, = T, (x), T* verifica 
4,1 e portanto é harmónica. 

Por outro lado, T* verifica as condições impostas 
emx=aey= b mas 


T*t<T emx=0 


a 


pelo que T* < T em todo o domínio. 
O produto de soluções dos problemas unidi- 


mensionais constitui pois um limite inferior. 
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4,1.3— A construção dum limite superior pode 
obter-se propondo 


0,693 É 
T=1-L(e É 4 
b À ; 


a qual verifica as condições aos limites em 
x=0€ey=0. 

Além disso T,< O (supondo a => b) em 
x=aey=b. 

Por outro lado 


dd g E me 


Pelo que podemos concluir que 
T<« 1, 
Combinando este resultado com o anterior : 


ú 0,693 x X 
a 


| y ! "a E yY/ | 
1-0)(1-S)<r<a- 2 (SE) 


) " 

A estimativa para T deve considerar-se ilustra- 
tiva. Tanto o limite superior como o inferior 
podem ser aproximados usando outras funções, 
cuja forma dependerá muitas vezes da intuição 
física que se tiver da solução. 


4.2 — Considere-se como segundo exemplo o 
caso de uma barra quadrada de lado 2 com fontes 
uniformemente distribuídas, Q=-=4, cuja fron- 
teira se mantém a T=0. 

Matematicamente será determinar T tal que 


4.2) vº T = para |x| <L |y|<| 


T=o0o para |x|=1ely/|/=1 


o problema é equivalente à determinação de u 
tal que 


4.3) xy CD 
u=x + y! x y CB 
Pois, como se verifica facilmente 


T=utx +y 
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Satisfaz 4.2) se u satisfizer 4.3, 
Admitamos como função de aproximação 


4,4) V = + ã| MV; E ã, Va 
onde 
v, = parte real de (x + iy)' 
v, == parte real de (x + iy)º 


sendo a, a, e a, constantes a determinar. 

As funções v, e v, são harmónicas, como o é 
a parte real e a imaginária de qualquer função 
analítica de variável complexa. 

Pelo teorema do máximo, o erro máximo e 
mínimo de tomar v por u verificar-se-á na fron- 
teira. O método consiste pois, simplesmente, 
em calcular |v-u| na fronteira e escolher 
a,, à;, a, de tal modo que Eu = | v-u | 
seja mínimo. 

Traçando a curva de v- u para cada lado (dada 
a simetria do problema basta neste caso traçar 
para um deles por exemplo, y==1) e conside- 
rando apenas v = a, +t- a, v, obtém-se, se 
a,=1,175ea,/=0,2 que | 1,75-0,2 v, (xº+1) | 
| <, 0,025 e portanto que | u-v | < 0,025 
Considerando agora v=a,+ta, v ta v € 
após ajustamento [5] se escolher v == 0,1786 — 
— 0,1801 v, + 0,006 v, verifica-se que em y =1 

| v-u | < 0,005 
que seria portanto o erro máximo cometido no 
domínio. A solução está portanto contida entre 


0,1786 — 0,1801 v, + 0,006 v, + x? + y? — 
— 0,005<. T<. 0,1786—0,1801 v, + 0,006 v, + 
+ x* + yº + 0,005 


Por exemplo, em x =y == 0 será 
1,1736 =. T (0,0) <. 1,1836 


4.3 — Os exemplos anteriores mostram suma 
riamente como os teoremas são susceptíveis de 
aplicação prática. O seu carácter elementar não 
revela porém completamente a sua potencialidade. 
Refira-se, em particular, a combinação com outros 
métodos de solução aproximada que se tenciona 
apresentar oportunamente. Estes facilitam a 
construção mais expedita de funções de 
aproximação da solução usada. Os teoremas 
de máximo permitem, através da fixação de 
alguns parâmetros dessas funções, a determinação 
dos limites do erro, aspecto basilar em qualquer 
solução aproximada. 
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TECNICA XVI 


é. D. U. 621.375.13:624,373.421 


REALIMENTAÇÃO EM SISTEMAS ELECTRÓNICOS 
LINEARES 


(Conclusão) 


4. CONFIGURAÇÕES TÍPICAS DE OSCILA- 
DORES DE REALIMENTAÇÃO 


4.1 Considerações gerais. Critério de Bar- 
khausen 


Quando, em relação ao esquema básico da 
fig. 2.1 estabelecemos a equação básica da rea- 
limentação 


A(s) 


“a 5 A 
(e) 1—b(s) A(s) 


(4.1) 


vimos ainda que, se se verificasse a condição, 


1=f(s). A(s) (4.2) 
o ganho total seria teoricamente infinito e o sis- 
tema seria portanto auto-excitado. Dissémos 
então que este era o caso dos chamados oscilado- 
res de realimentação. 

Efectivamente a condição (4.2) corresponde 
aos zeros da função complexa de variável com- 
plexa 

E(s)=1—F(s). A(s). (4,3) 

De uma maneira genérica, um oscilador de rea- 

limentação pode ser representado pelo esquema 


AMPLIFICADOR 
(ELEMENTO ACTIVO) 


MALHA 
PASSIVA, DE 
REALIMENTAÇÃO 


Fig 4.1 


básico da fig. 4.1 que corresponde ao esquema 
da fig. 2.1 no caso de ser y(t)=0ec(t) =r (t). 
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pelo 


Prof. ANTÔNIO DE CARVALHO FERNANDES 


A equação (4.2) equivale à condição (4.4) que 


io sá 
= 


(4,4) 
corresponde ao conhecido Critério de Barkhausen (*). 
Deve notar-se que o significado de G(s) e A(s) 
continua a ser o mesmo das equações (2.1) e (2.2) 
tendo em vista que o ganho A (s) será, em geral, 
afectado pela presença da malha de realimen- 
tação. 

Podemos pois dizer que há dois elementos 
essenciais num oscilador de realimentação: a) 
um elemento activo que, para pequenos sinais, 
terá comportamento linear mas que, para sinais 
de uma certa amplitude apresentará uma acção 
de não-linearidade que tende a limitar essa am- 
plitude; b) uma malha de realimentação consti- 
tuída por elementos passivos. Quando estes 
elementos são combinados num sistema de reali- 
mentação positiva — a que corresponde a condição 

F(s)|=/1—B(s). A(s)|<1 (4.5) 
o sistema tende a ser auto-excitado e qualquer 
regime livre acidental, que tenha lugar no circuito, 
levará a oscilações de amplitude crescente. 
Quando a não-linearidade do elemento activo 
influenciar o seu comportamento de tal modo 
que as oscilações se mantêm com amplitude cons- 
tante, os parâmetros do elemento activo e da 
malha de realimentação serão tais que, para uma 
certa frequência, se verificará a equação (4.4). 

É de notar que, no início das oscilações, o ele- 
mento activo tem funcionamento essencialmente 
linear e que é a sua não-linearidade que limita a 
amplitude das oscilações. Pode portanto pergun- 
tar-se; e em regime permanente, o funcionamento 
é ou não linear? De facto, os dispositivos dos 


(*) O critério de Barkhausen foi estabelecido para 
o caso de osciladores usando válvulas electrónicas e em 
condições mais restritas. 


175 


osciladores electrónicos podem ter várias classes 
de funcionamento: em classe À, a operação será 
linear e nas classes B ou C o funcionamento dos 
dispositivos electrónicos será essencialmente não- 
-linear. No nosso estudo, vamos, contudo, admitir 
que o funcionamento é linear e é nesta hipótese 
que é válida a equação (4.4). 


4.2 Configurações típicas 


Às quatro configurações típicas de amplifica- 
dores com realimentação corresponderão, respec- 
tivamente, quatro configurações idênticas para 
osciladores de realimentação ; trata-se essencial- 
mente das quatro maneiras básicas de associar 
o Amplificador (ou elemento activo) e a Malha 
de Realimentação. Os respectivos esquemas 
básicos são os seguintes (figs. 4.2 a 4.5), em que 
supomos um regime sinusoidal e admitimos 
portanto que as tensões e as correntes são funções 
da variável s = jm. O Cálculo Operacional dá 
assim lugar ao Cálculo Simbólico em que as 
tensões e as correntes são representadas pelas 
respectivas amplitudes complexas. 

A) Ligação série (na saída) — série (na entrada) 


ELEMENTO 7 
ACTIVO 


MALHA 
PASSIVA 


Fig. 4.2 


B) Ligação derivação (na saída) — série (na 
entrada) 
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C) Ligação derivação (na saída) — derivação 
(na entrada) 


ELEMENTO 
ACTIVO 


MALHA 
PASSIVA 


Fig. 4.4 


D) Ligação série (na saída) — derivação (na 
entrada) 


ELEMENTO 
ACTIVO 


MALHA 
PASSIVA 


Fig. 4.5 


4.3 Critério de 4 = 0 


Em qualquer destas configurações (figs. 4.2 a 4.5) 
não temos excitação exterior e, portanto, as 
equações relativas ao sistema são necessáriamente 
equações lineares e homogéneas. Para que este 
sistema tenha soluções não nulas, é preciso que 
o respectivo determinante seja nulo. 

à =0 (4.6) 

Este é um dos critérios mais simples para a deter- 
minação das condições de oscilação. Atendendo 
a que os parâmetros do circuito serão números 
complexos (funções de jw), a condição (4.6) 
traduz-se assim, em geral, por duas equações. 

Se, nos esquemas das figs. 4.2 a 4.5, quisermos 
admitir que se trata de osciladores sem carga 
exterior, teremos de considerar: Z, = nas 
configurações em que a ligação na saída se faz 
em série—caso das figs. Ze 5;e Z, = nas 
configurações em que a ligação na saída se faz 
em derivação — caso das figs. 4.3 e 4,4, 
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44 Exemplo 


Como exemplo da ligação derivação-derivação 
(fig. 4.4), sem carga exterior, consideremos o 
esquema simplificado da fig. 4.6 em que o 
elemento activo é um andar de amplificação, com 
um transistor numa montagem de emissor comum, 
e a malha passiva de realimentação é formada 
pelas impedâncias Z,, Z, e Z. que supomos serem 
reactivas puras. 


| o — À —— —o——e— 
| | | 
GE CO CR NE 

Fig. 4.6 


Podemos ainda notar que, em relação ao 
esquema da fig. 4.4, a fig. 4.6 poderia tomar o 
aspecto representado na fig. 4.6-a o que equivale 
a admitir que Z, =Z, e que a malha passiva de 
realimentação é formada apenas pelas impedân- 
cias Z, e Z,. 


O— — =. 


No Apêndice N.º 5 fazemos o estudo deste 
circuito e determinamos as condições de oscilação 
(em regime permanente) quer a partir do critério 
do determinante nulo — equação (4.6) — quer a 
partir do critério de Barkhausen — equação (4.4). 
Admitindo que 


Z=]X B=jX, Z=]Xs (4.7) 
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conclui-se que deverá ser 


Xj + X + X,==0 (4.8) 


X = he X (4.9) 
Isto significa que as impedâncias X, e X, serão 
da mesma natureza, por exemplo, capacitivas 
(como no oscilador Colpitts) e, nesse caso X, 
deverá ser indutiva. O esquema simplificado do 
Oscilador Colpitts utilizando uma montagem de 
emissor comum é dado na fig. 4.7 em que 


| 1 
1 
Z,=)X,=—) 1 C, (4.11) 
Z,=jX,= jeL (4.12) 
] 
| =) Ea E 
l; 
Fig. 4.7 
devendo ser, portanto, 
a. Q+C 
L C, E, (4.14) 


Estas condições definem-nos a relação entre as 
capacidades, para que o sistema possa oscilar, 
e dão-nos a frequência angular de oscilação; é 
de notar que esta é a frequência própria do cir- 
cuito oscilante do esquema da fig. 4.7, que supu- 
semos sem perdas. 

Muitos outros exemplos poderiam ser citados 
tendo em consideração as diferentes configura- 
ções possíveis para o andar de amplificação (por 
exemplo de emissor comum, de base comum e de 
colector comum, no caso de transistores), para a 
malha passiva de realimentação (onde é possível 
uma vasta gama de configurações e para a ligação 
entre estes dois elementos (como ilustrado nas 
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figs. 4.2 a 4,5). Nem sempre o conjunto escolhido 
destas configurações é compatível no sentido de 
poder ser auto-excitado e constituir assim um 
oscilador. 


5. ESTABILIDADE DE SISTEMAS LINEARES 
COM REALIMENTAÇÃO 


5.1 Definição de estabilidade 


De um modo geral, um sistema estável perma- 
necerá em repouso excepto quando estiver sobre 
a acção duma excitação exterior e voltará à 
situação de repouso se essa excitação for retirada. 
Matemáticamente, a estabilidade pode ser definida 
com precisão a partir da resposta à função im- 
pulsiva unitária: — um sistema diz-se estável se a sua 
resposta à funcão impulsiva unitária tender para zero 
quando a variável tempo tender para o infinito. 


5.2 Resposta à função impulsiva unitária 


Designamos por k (t) a resposta, no domínio 
da variável t, à função impulsiva unitária 5, (t) 
cuja transformada de Laplace é igual à unidade. 
Neste caso, verifica-se que k (t) é a transformada 
inversa de Laplace da função de transferência 
K (s) do sistema, como se mostra no Apêndice VI. 

O problema da determinação da estabilidade 
de um sistema reduz-se pois ao problema de 
averiguar se a respectiva função k (t) (resposta à 
função impulsiva unitária) tende para zero 
quando t tender para infinito. Se considerarmos 
um sistema com realimentação em que a função 
de transferência global (ou ganho global), dada 
pelas equações 


X(s) Als). 


K (s) = = = 
Y(s) 1-É(s) A(s) 


RE (1) 


o problema resumir-se-ia em determinar a res- 
pectiva tranformada inversa k(t) e averiguar, 
directamente, se tende para zero quando t tender 
para infinito. 

Este procedimento implicaria, contudo, a deter- 
minação da transformada inversa de K(s),o 
que pode ser difícil e laborioso. 
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5.3 Análise da função de transferência. Cri- 
tério de Routh 


Por isso, têm sido propostos vários critérios 
com vista à determinação da estabilidade de um 
sistema pela análise da sua função de transfe- 
rência global. Esses critérios têm, necessáriamente, 
em consideração as propriedades da transformação 
inversa de Laplace. 

No Apêndice VII mostramos que a transfor- 
mada inversa de uma função K (s) terá termos 
afectados do factor e 7! sendo 7 => 0 quando a 
função tiver polos (simples ou múltiplos, reais ou 
complexos) situados na metade da direita do 
plano da variável complexa s = + jm. Isto 
significa que, nessas condições, o funcionamento 
do sistema não é estável. 

Obtemos assim a conclusão fundamental de 
que o sistema será instável se existirem polos da 
função de transferência global dada pela equação 
(1) ma metade da direita do plano da variável 
complexa s==7 + jw, Pode ainda provar-se que 
o sistema será estável se, e só se, não existirem, 
polos de K (s) na metade da direita do plano da 
variável s. O mesmo é dizer que o sistema 
será estável se, e só se, não existirem zeros da 
função [1 — É (s) A(s)] na metade da direita do 
plano de s. Esta conclusão constitui o chamado Cri- 
tério de ROUTII no domínio do Cálculo Operacional. 

Daqui se conclui uma regra prática: deter- 
minar os zeros de [1 — É (s) A(s)| e verificar se 
satisfazem a condição referida. 

Nem sempre, contudo, é fácil fazer esta deter- 
minação. Por outro lado, em circuitos práticos, 
nem sempre é fácil obter as expressões analíticas, 
de C(s) e A(s) mas é algumas vezes possível 
determinar os seus valores (em módulo e argu- 
mento) experimentalmente, em regime sinu- 
soidal (s==jm), para diferentes valores da 
frequência. Daqui, a ideia de, a partir do diagrama 
da função [1—L(jw)A(jw)] para valores de 
w variando de — coa + co, tirar conclusões quanto 
a estabilidade do sistema. Vejamos que o chamado 
Critério de NYQUIST nos dá satisfação a este respeito. 


5.4 Critério de Nyquist 
O criterio de NYQuIsT, original, aplica-se ao 
caso especial de não haver zeros da função K (s) 


(ou polos de [1—fL(s)A(s)|) na metade da 
direita do plano da variável complexa. Este caso 
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especial abrange, contudo, a maior parte dos 
amplificadores com realimentação (para os quais 
o critério era originalmente destinado) e também 
a grande maioria dos sistemas com realimentação 
que se utilizam na prática. Em relação à função 
de transferência global, dada pela equação da 
realimentação (1), o critério de NYQUIST pode ser 
enunciado como segue: 

Admitindo que não há polos de [1 — É (s) A (5) |] 
na metade direita do plano de s— o que se veri- 
fica na quase generalidade dos sistemas práticos 
com realimentação — uma condição necessária 
e suficiente, para a estabilidade do sistema, é que 
o contorno fechado que constitui o lugar 
geométrico, no plano (u--jv) da função 
[1—B(jw)A (jv)], tomado sobre toda a gama 
de frequências positivas e negativas, não envolva 
a origem dos eixos. 

Um enunciado equivalente será ainda: 

O lugar geométrico da função [É (jw) A (jw)] 
não deve envolver o ponto (1,0), no plano 
(u -jv), quando w for tomado sobre todos os 
valores positivos e negativos. Este segundo enun- 
ciado corresponde, evidentemente, ao primeiro, 
se considerarmos uma mudança na posição da 
origem dos eixos. 

No Apêndice X dá-se um exemplo prático 
de aplicação deste Critério. 


A demonstração do critério de NYQUIST envolve 
os seguintes passos essenciais: 


a) A condição de estabilidade obriga, como já 
vimos, à função [1 — É(s) A (s)] não ter zeros 
na metade da direita do plano da variável s. 

h) Se uma curva fechada, no plano de s, envolver 
um número Z de zeros e um número P de 
polos de uma função F (s), o contorno fechado 
que lhe corresponde no plano da função 
F(s) = u--jv envolverá a origem dos eixos 
tantas vezes quanto for o valor diferença 
Z—P. Esta afirmação demonstra-se a partir 
de um teorema clássico da teoria das funções 
de variável complexa, como se mostra no 
Apêndice VIII. 

c) O contorno, no plano de s, limitado pelo 
eixo imaginário e por uma semicircunferência 
de centro na origem dos eixos e de raio 
infinito (chamado o Caminho Fechado de 
NYQUIST) abraçando portanto todos os polos 
e zeros de uma função F(s) na metade da 


TECNICA Nº “91 


direita do plano de s, transforma-se, no 
plano de F(s) == u + jv no lugar geométrico 
da função F(jm) isto é, de F(s) quando se 
toma s =— |», considerando todos os valores 
positivos e negativos da frequência. Esse 
lugar geométrico designa-se por Diagrama de 
xyquisr. No Apêndice IX faz-se a demons- 
tração desta afirmação. 

d)Se a função referida em bh) e c) for 
[1—B(s) A(s)] e se o contorno no plano 
u-jv que representa o lugar geométrico de 
[1-E(jwo) A(jo)] para—-co<So<<+co 
não envolver a origem dos eixos, isso quer 
dizer que Z — P == 0. Como admitimos que 
P-==0 (na quase generalidade dos sistemas 
com realimentação), conclui-se que será 
também Z=-0 e que não haverá zeros 
daquela função na metade da direita do plano 
de s e, portanto, que o sistema é estável. 
Se o contorno envolver a origem, uma ou 
mais vezes, então será: 


7 sp (2) 
e, por ser P = 0, será 
Zoo 1 (3) 


o que significa que haverá, pelo menos, um 
zero na metade da direita do plano de se o 
sistema será instável. 

Esta é pois a justificação do Critério de NYQUIST 
a que nos referimos. 

É de notar que o lugar geométrico dos pontos 
da função [1—LB(jw) A(jm)] no plano 
u--jv se pode obter experimentalmente 
a partir de medidas de É e A (os seus 
módulos e os seus argumentos), sobre uma 
larga faixa de frequências; bastará teorica- 
mente considerar a faixa (0, + o2) pois que 
haverá simetria conjugada neste diagrama, 
o que quer dizer que o gráfico relativo a 
—co<w< 0 é a imagem de espelho, em 
relação ao eixo horizontal, do gráfico relativo 


aan os, 
5.5 Outros critérios 


Os critérios referidos de ROUTH e NYQUIST, 
são muitas vezes apresentados de forma diferente 
da que vimos. Têm sido ainda apresentados vários 
outros critérios para a determinação das condições 
de estabilidade. 
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APÊNDICE V 


CONDIÇÕES DE OSCILAÇÃO. OSCILADOR 

DE REALIMENTAÇÃO EM LIGAÇÃO DERI- 

VAÇÃO-DERIVAÇÃO. CRITÉRIOS DE DETER- 
MINANTE NULO E DE BARKHAUSEN., 


Em relação ao esquema da fig. 4.6 (do Cap. 4), 
podemos traçar o esquema equivalente da fig. VI, 
numa versão muito simplificada do esquema do 
transistor, onde 1, e 1, representam amplitudes 
complexas, em Cálculo Simbólico, e são funções 
da frequência angular ». Também as impedâncias 
serão função de w, As equações relativas ao 
circuito são: 


Lhe +(L+1)Z,=0 (1) 


ds 


(1, + 1,)Z,+ 1,2,;— (ho 1—1,)Z2,=0 (2) 
a que corresponde o sistema: 


1 (he HZ) 
LtZ,;—h5 25) 


à LT; = 
+1L,(Z,+Zo + Z;))=0 


(3) 


A condição de determinante nulo é portanto: 
(Me AU Mp L+L)-ZUIMZ—hZ,)=0 
(4) 
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Admitindo que 
L=jÃX, Lo=jX, Z;=jX; (5) 
obteremos 


BilA, + X + AÃ TI [X, (X, + X + X5)— 
=X, +h,X,X,|=0 (6) 


Esta equação será satisfeita quando se verifi- 
carem as condições 


X + XD (7) 
X,=hçX; (8) 
que representam assim, as condições de oscilação. 


Em relação ao esquema da fig. 4,6-a do Cap. 4, 
o respectivo esquema equivalente será: 


b= E Ela = . nes 
Ne |] a | 
' | | 1 | | 1 
O * .— O + em 
Fig. V.2 
Vi 
Designando por A = T (9) 
p=, (10) 
Vo 
o criterio de Barkhausen, generalizado, será 
1 
E=—.. (11) 
A 
Em relação a este esquema será 
V=I he (12) 
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U=—h(he+14 e jz, (13) 
A=— É (he 414 E (14) 
Nie l 
h;. || Z 
2 EE a e Va (15) 
hj |Z,+2Z; 
E hj | Z, qa he Zi À 
he lZ, TZ; heZytheZ;t2Z,Z, 
(16) 


e 1 pass : 
A condição É = E leva assim à equação 


bu li, E) SALA AZ 
(17) 
Fazendo 


Zi= já; Lo =X Zs= 5X; (18) 


sera 


hi !X, + X, + X)—XX(1+h)—X X=0 


(19) 
que será satisfeita quando for 
X+X,+X;=0 (20) 
X(1-+Hh)=—X, (21) 
ou ainda quando for 
X,+X,+X,=0 (22) 


X,=hçX, 


como já tínhamos encontrado a partir do critério 
de determinante nulo. 


APÊNDICE VI 


RESFOSTA DE UM SISTEMA, À FUNÇÃO 

IMPULSIVA UNITÁRIA. TRANSFORMADA 

INVERSA DA FUNÇÃO DE TRANSFERÊNCIA 
DO SISTEMA 


Designemos por k (t) a resposta, no domínio 


da variável t, à função impulsiva unitária e 
verifiquemos que k (t) é a transformada inversa 
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de Laplace da função de transferência K (s) do 
sistema, 

Para concretizarmos melhor o nosso estudo, 
imaginemos um sistema físico cujo comporta- 
mento é traduzido pela seguinte equação dife- 
rencial linear 


| i 
Ba SU =y() (1) 


== 0 


em que y (t) é função excitação aplicada ao sis- 
tema e x(t) é função resposta; os coeficientes 
a; são constantes conhecidas. Se apliçarmos a 
transformação de Laplace a esta equação obtere- 
mos, no domínio da variável complexa a equação 


X()=Y(s5) K (s) + Q(s) K (5) (2) 


em que K (s) é a chamada função de transferência 
do sistema e Q (s) é uma função que depende 
das condições iniciais. A função de transferência 
é dada pela expressão 


K()=—>—— 


2 a, s' (3) 


No caso particular de ser n==2, os valores 
de K (s) e OQ (s) serão, respectivamente 


1 


K (9) = ——— 
atasta,s 


(4) 
O (s)=(a,s+ta,)x(0) + a, x" (0) (5) 


em que x(0) e x (0) são respectivamente, os 
valores de x (t) e da sua derivada x' (t) no ins- 
tante inicial t = 0. 

Se admitirmos como condições iniciais x (0,=0, 
x (0) =0 .. gt (0) = 0, isto é, se admi- 
tirmos que o sistema está inicialmente em repouso 
e se admitirmos ainda que há continuidade na 
origem da variável e das suas derivadas até à 
ordem (n-1), será Q (s) = 0 e a função de trans- 
ferência K (s) corresponderá ao cociente 


X (5) 
K == — 
O ro) 


(6) 


como temos considerado no nosso estudo. 


Analisemos agora 0 caso particular de a função 
de excitação y (t) ser a função impulsiva unitária 
na origem 9, (t) cuja transformada de Laplace é 
igual a 1 


Ly(0M)=2f, (0]=1 (7) 
Neste caso será 
K (s)=Xs (s) para y(t)= 0, (t) (8) 


No domínio do tempo, se designarmos por 
k (t) a transformada inversa da função de trans- 
ferência, 


k (t) =" AK (9)) (9) 
concluiremos que 
k (t) = x (t) (10) 


Isto significa portanto, que a resposta de um sis- 
tema nas condições referiãas, para uma excitação soh a 
forma da função impulsiva unitária a (t), é igual à 
transformada inversa da função de transferência do 
sistema. Por isso k (t) é designada por resposta à 
função impulsiva unitária. 


APÊNDICE VII 
LOCALIZAÇÃO DOS POLOS DA FUNÇÃO 
DE TRANSFERÊNCIA GLOBAL, COMO CRI 

TÉRIO DE ESTABILIDADE DO SISTEMA 


Consideremos a função de transferência global, 
dada pela equação básica da realimentação : 


Et) = — SE (1) 


A respectiva transformada inversa será, pelo 
teorema da inversão 


1 lim f 
67 Pes 


Y=ib 


ar 
E | 15 


k (t) = K (s)el ds (2) 


em que K (s) deve ser uma função complexa, de 
variavel complexa, analítica no semiplano 


Re [s| > 7%. (3) 
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O valor de y, no integral de inversão deve 
pois ser escolhido de modo que todas as singula- 
ridades de K (s) existam no semiplano Ra [8] <<? 

Admitamos agora que a função K (s) tem 
como polos os pontos s,,s,...s, no plano da 
variável complexas. Estes serão também os polos 
da função K (s) e”. Consideremos ainda que os 
residuos da função integranda, no integral de 
inversão, K (s) e*!, são, respectivamente R (s,), 
Rio) cota) 


Pode provar-se que, sendo c uma curva fechada 
formada por um segmento de recta 9 == com 
os extremos nos pontos y + jb e y—jfb e por 
um arco de curva qualquer, de tal modo que no 
interior da região limitada por c existam todos 
os polos de K (s), se verifica a equação 


Trib i 
lim K (5) Rd ds ——— ) K (5) e“ ds . (4) 
rias Y—iê A 


Tendo agora em consideração o Teorema dos 
Resíduos de Cauchy 


| F(s) ds=27; > R(s) (5) 
c j=l 


podemos escrever que a transformada inversa 
de K (s) será 


m 


kt)= 5 R(s) (6) 
| =] 
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que traduz o Método dos Resíduos para a deter- 
minação da transformada inversa de Laplace. 
Este método consiste em determinar os resíduos 
da função K s) e* e obter o somatório como se 
estabelece na equação (6). Convém ainda referir 
que, se uma função F (s) for dada pelo cociente 
de duas funções analíticas 


| P (9). 
E (s) = (7) 
Q (s) | 
o resíduo num polo simples s, é igual a 
P (s,) 
R (89) = —— (8) 
| OQ (s,) 


e que, dum modo geral, o resíduo num polo 
múltiplo, s, de ordem N, é igual a 
| 8 = 5 q 


1 aNo' 
(9) 


(N—1)! | gsN- 

Temos assim uma regra prática para a deter- 
minação dos resíduos da função F(s) =K (s) e! 
e da transformada inversa k (t). | 

A determinação dos resíduos da função F (s) 
é facilitada quando a função K (s) é dada pelo 
cociente de dois polinómios 


R(s)= (s-s)N E(9) 


Et) — 289. (10) 


em que p(s) e q (s) são polinómios sem factores 
em comum, sendo o grau de p (s) inferior ao de 
q (s). Os polos de K (s) serão agora as raizes de 
q (s). A determinação dos resíduos que consti- 
tuem as parcelas do somatório que nos dá a 
transformada inversa levam às seguintes regras 
de aplicação prática de que tiraremos conclusões 
quanto ao problema em jogo de a função k (t) 
dever tender para zero quando t tender para 
infinito, para que o sistema seja estável. 

a) Por cada factor (s-s,)), não repetido, no 
denominador, isto é, para cada polo simples s, da 
função K (s) haverá o termo 


8 (11) 
q Is) 


na transformada inversa k (t). 
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Para valores de S, reais, conclui-se imediata- 
mente que deverá ser 


s5<O (s; real) (12) 
para que a resposta à função impulsiva unitária 
seja atenuada exponencialmente no tempo. 

Para valores de 5,9, + jw, complexos, conclui- 
-se também que 

3, = Re [8;] <0 (13) 
Em conclusão, para que o sistema seja estável, 
a função K(s) deverá ter todos os seus polos simples na 
metade da esquerda do plano da variável complexa defi- 
nida por Re [s]< 0. 

Se nos debruçarmos sobre o caso de K(s) ter 
dois polos imaginários simples + jw, e — jm, sen- 
do, portanto, 7,==0, haverá dois termos em que 
nos aparecem os factores e jojte emiot , que 
correspondem a um regime sinusoidal de fre- 
quência angular »,. O sistema não é estável mas, 
neste caso, diz-se, algumas vezes, que o sistema é 
marginalmente estável. 

Se, além dos factores el“ite e nos 
aparecesse um factor e “if em que > O isso 
significaria que as oscilações tenderiam a aumentar 
de amplitude e o sistema seria instável; sendo 
5,< O, as oscilações tenderiam a diminuir de 
amplitude e o sistema seria estável. 

b) Por cada factor (s — s,)" no denominador de 
(7), isto é, por cada polo de ordem r da função 
K (s) haverá os seguintes r teremos em k(t) 


— jwt 


r "e da 
A Ca (14) 
pe (2=L)! (11)! 
em que 
di(s)=K(s). (s—s,)" (15) 


Daqui se conclui que o factor e *;* torna gené- 
rica a regra de a função K'(s) só dever ter polos, quer 
simples, quer múltiplos, na metade da esquerda do plano 
da variável complexa s, para que o sistema seja estável. 

A conclusão referida é válida quer os polos 
s; correspondam a números reais, quer a núme- 
ros complexos. 
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APÊNDICE VIII 


CONTORNO FECHADO NO PLANO DE 5, 

ENVOLVENDO z ZEROS E p POLOS DE F (s). 

CONTORNO CORRESPONDENTE NO PLA- 

NO DE F(s) ENVOLVENDO (z— p) VEZES 
A ORIGEM DOS EIXOS 


Consideremos o seguinte teorema clássico : 


Seja F(s) uma função analítica e uniforme num 
domínio R; seja C uma curva fechada simples, 
em R, tal que F(s) não é igual a zero em qual- 
quer ponto sobre o contorno C e o interior de € 
contém, um número finito de singularidades deF(s), 
todas elas polos. Se o número de zeros no 
interior de C (contando um zero de multiplici- 
dade m, m vezes) é z, e se o número de polos 
de F(s) no interior de €C (contando um polo de 
ordem n, n vezes) é p, então 


T TED a 
— pipe (1) 


Se reparamos que 


; x ” 
dlog F(s) e: (s) (2) 
ds Fís) 
e se escrevermos F(s) sob a forma 
F(s) = | F(s) | e!"t E (3) 


LU 
sera 


logF(s)=log | F(s) | + jarg F(s) (4) 


APÊNDICE 1X 


CAMINHO FECHADO DE NYQUIST (NO 
PLANO DE s) DIAGRAMA DE NYQUIST 
(NO PLANO DE F (s)) 


O caminho fechado de NYQUIST é um contorno 


fechado, no plano da variável complexa s = q + jw, 
que abrange completamente a metade da direita 
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Podemos assim escrever que, de acordo com o 
teorema referido, será 


1 Ff rig). 1 | 
p= 2 (EO faftog Fo) = 
ii Di] o i sr | [og F(s)| 

o Cc 


= A [og E(s)] - 


2) 


am E Host | 4 je ted] | e 
2 T) C 
1 

= o [og [E (9) 1 + 
27] C 


1 
E 
fa 91, (5) 


A primeira destas duas parcelas é nula porque 
log | F(s) | volta ao seu valor original quando 
circulamos uma vez ao longo do contorno 
fechado C. 


Será pois 


z-p=—E [arg E(s)] (6) 
2 C 


Como já vimos, ao contorno fechado C, no 
plano de s, corresponderá um contorno tam- 
bém fechado, relativo a F(s) no plano de 
F(s)=u + jv, de tal modo que a variação no 
ângulo do arg F(s) é (2z— p) vezes 27, isto é, tantas 
vezes 27 quanto é o valor n==z—p. Isto quer 
dizer que o contorno no plano de F(s) envolverá 
a origem dos eixos uma vez se forn=1 en 
vezes se for n==n, Se o contorno de F(s) não 
envolver a origem, será n=-0. Em resumo, 
o contorno de F(s) rodeia a origem dos eixos 
n==(z—p) vezes, como queríamos provar. 


do plano de s. É formado, da maneira que se 
ilustra na fig. IX.1, em que o raio R do semi- 
círculo da direita tende para infinito e o raio ? 
dos pequenos semicírculos, à volta dos polos de 
F (s), situados sobre o eixo dos ww», tende para 
zero. Esses semicirculos são utilizados para evitar 
que o contorno fechado contenha singularidades 
da função F (s). Este caminho fechado abrangerá 
assim todos os polos (em número p) e todos os 
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zeros (em número z) da função, situados na 
metade da direita do plano da variável s. 


- ]Wp) 


* Folos possíveis da 
função FIs) | 


Fig. IX.1 


Para obtermos o contorno correspondente ao 
caminho fechado de yquisr, no plano de 
F(s) =u-tjv precisamos considerar três tipos 
de troços no caminho fechado de NYQUIST: a) as 
semicircunferências de raio infinitesimal à volta 
de polos situados no eixo dos w»; b) a semi- 
circunferência infinita; c) os troços situados no 
eixo dos wo, 

Consideremos em primeiro lugar que existe 
um polo na origem dos eixos no plano de s, 
isto é, para s == O como se ilustra na fig. IX.2. 
Os pontos situados sobre essa circunferência 
podem ser expressos por 


s, = dr el! (1) 


Fig. IX? 
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Na vizinhança de um polo s,, o valor de F (s) 
é determinado, essencialmente, pelo valor do 
polo e para um polo múltiplo de grau N, a 
função F(s) na sua vizinhança pode escrever-se 
dentro de certa aproximação 


1 
(s—s,)º 


F (s) = (2) 


em que 1 é constante. 
Se o polo estiver no eixo dos mw será s, = j9,, 
e poderá escrever-se 


F (s) = 


(s— jo) 


(3) 


Se o polo estiver na origem dos eixos será 
mw ==0 e os valores de F (s) correspondentes à 
semicircunferência de raio infinitesimal serão 


F(s)=1 dr" EN (4) 


Lj 7 
a— — . O respectivo 
2 2 


em que x varia de + 


lugar geométrico será um contorno circular de raio 
infinito correspondente a um ângulo igual a Nr 
quando w varia desde um valor positivo muito 
pequeno w = O” até um valor negativo muito 
pequeno v= 0. 

Consideremos o caso de dois polos simples 
S=j” e ss=-—j»,a que corresponderia, no 
denominador de F(s) o factor (sº + «,2). 

Suponhamos, para maior simplificação, que 


1 1 
fimo É ul É fá 
2 + E) 2 (s + j 60,) (s —— ] nm) 


O lugar geométrico dos pontos relativos a F(s) 
no plano de F (s)=u + jv correspondente ao cami- 
nho fechado na fig. IX.3 é representado na fig. IX.4, 
em que os troços correspondentes às semicircun- 
ferências à volta de w», (desde ww,” a w” e desde 


— 6 + a —tw corrrespondem a semicircunfe- 


Fig. IX.3 


rências de raio co. Pode pois concluir-se que as 
semicircunferências de raio infinitesimal, que 
semi-rodeiam os polos no eixo dos w» (no plano 
de s), dão origem, no plano de F (s) a percursos 
circulares de raio infinito. 

Consideremos agora a semicircunferência de 
raio infinito que faz parte do caminho fechado 
de xyQuisrT. Para valores de s sobre essa circun- 
ferência, os valores de F(s) poderão ser expressos, 
dentro de certa aproximação por 


dd Cor uted  i (6) 


em que z é o número de zeros da função F (s) 
e p é o número de polos. 

Nos sistemas físicos correntes, p excede z. 
Sobre a semicircunferência referida será 


s=R.e (7) 


em que R>0O e O variade—— a + — 
à 2 
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(? à 
, N 
f ds ' 
j | 1=1) ] 
E E DE ms ==, 
RR: á ' Es (- E - 
P O , = 
dn w=+e | ! cos 
! 1 f 
N | wi | st 
A La =| = 
“ | Ed 
o 
| 
Vig. IXA 
O lugar geométrico de (8) 
z-p  Azp)0 (8) 


F(s) =1.R e) 


que lhe corresponde, no plano de F(s) será de 
forma circular e raio infinitesimal, a que corres- 
ponderá, no limite, um ponto que será a origem 
dos eixos. 

Vemos, portanto, que no caminho de NYQUIST, 
as semicircunferências de raio infinitesimal dão 
origem, no contorno correspondente de F's), a 
percursos circulares de raio infinito; e que a 
semicircunferência de raio infinito dá lugar a um 
ponto que estã na origem dos eixos. 

O lugar geométrico dos pontos de F (s) corres- 
pondentes aos pontos no eixo dos w», no plano 
de s, isto é, para 


s= 0 


<Lm<os (9) 


será o chamado diagrama polar de F(jm»). Este 
diagrama polar apresenta simetria conjugada, o 
que quer dizer que o gráfico relativoa— co. m<.0 
é a imagem de espelho, em relação ao eixo hori- 
zontal, do gráfico relativo a O <w < 2, Trata-se 
de um contorno fechado e é a este lugar geomé- 
trico que se dá a designação de Diagrama de Nyquist. 
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APÊNDICE X (*) 


EXEMPLO DE APLICAÇÃO DO CRITÉRIO 
DE NYQUIST 


Como exemplo de aplicação do critério de 
Nyquist vamos estudar o caso simples do mul- 
tivibrador astável, empregando válvulas elec- 
trônicas (**). 

O multivibrador astável pode ser considerado 
como um amplificador de dois andares de acopla- 
mento RC em que a saida se encontra directa- 
mente ligada à entrada, proporcionando assim 
uma realimentação positiva de 100'/,. Um esquema 
prático de multivibrador encontra-se répresen- 
tadonafig. X.1. Para maior simplicidade, supõe-se 
em relação a este esquema que os dois andares 
associados em cadeia são idênticos. 

Para avaliar a estabilidade deste sistema recor- 
rendo ao critério de Nyquist, bastará de acordo 
com o enunciado dado em (5.4) verificar se o 
lugar geométrico da função A(jw). B(jm) 
engloba o ponto de coordenadas (1,0) no plano 
u--jv quando w for tomada sobre todos os 
valores positivos e negativos. 


Fig. X.1 
Vera 
V 


é, neste caso, independente da frequência e igual 

a + 1. Não teremos então mais do que traçar o 

diagrama polar de A(jw) no plano (u+jv). 
Como os dois andares são idênticos, 


Da fig. X.1 podemos concluir que £ — 


o 


A(jo)=A,. A, =(Aj) (1) 


(*) Este Apêndice foi preparado pelo Assistente Eng.º 
Pedro Luís Borges Teixeira. 

(**) Este circuito encontra-se estudado respectivamente 
para o caso de válvulas e transistores nos 5 10.16 do 
Cap: X e 17.8.7 do Cap. XVII do CURSO DE ELECTRÓ- 
NICA APLICADA, do Autor (Ed. AEIST, 1967). 
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Para a determinação aproximada da variação 
do ganho de tensão com a frequência em relação 
a um dos dois andares, consideremos o esquema 
equivalente da fig. X.2 


Em relação a este esquema, podemos, de acordo 
com a equação (13) do Apêndice II, escrever 
para expressão do ganho de tensão A, na faixa 
central das frequências 


Va 
A, = v. = — gm Ro (2) 
Reg =", //R$//R (3) 


Na zona inferior da faixa, não poderemos 
desprezar a reactância do condensador de acopla- 
mento C,, e análogamente à equação (15) do 
Apêndice II podemos escrever : 


1 
A, = À, ==——— (4) 
1—-j 
ta) 
com 
| = (5) 
RIC | 
e 
ro Rj 
R=R4+—— (6) 
* +R 


Para a zona superior da faixa de frequências 
podemos desprezar a reactância de C,, mas não 
poderemos deixar de considerar o efeito das 
capacidades interelectródicas da válvula do 
andar 1, assim como o da capacidade de entrada 
do andar 2. À semelhança da equação (23) do 
Apêndice II, teremos para expressão do ganho: 


RE (7) 


187 


com 


1 | 
tj E E (8) 
Es Ra 


C. =—— Es ), + C,, (9) 
Sendo Req dada pela expressão (3). 
A amplificação de tensão global para as fre- 


quências inferiores da faixa será por conseguinte 
dada de (4), (2) e (1). 


“ —8Bm Ran : Bu Rar 12% 
ara] = a PT omite (10) 
ad les 1+( ) 


" 


E para as frequências superiores: 


“— Em Rea E Bm l Rea —124 
(A)s — | = | / tú | e 
Ea a ) 
iria 1 + 0) (11 
bh] 
com 7 == arctg — (12) 
(a) 
” Eri 
e 9" == arctg — (13) 
Do 


A partir das equações (10), (11), (12) e (13) 
podemos traçar, a partir do quadro X.3 a variação 
aproximada de A (e por conseguinte de AÉ) com 
a frequência. 


| o | ul wa | DO 
mx [E] =To 
| E 
2% — | — Ea | te 
| Ea | Z | 
| | E E— 
| | 
eg o (EB)  |o= 
k o DT | 
| | FEET Nas 
(A). a = em “eq 0 
= 


Quadro X,3 


A fig. X.4 apresenta, de acordo com o quadro 
X.3 o diagrama polar de Aú no plano u + jv. 
Para a faixa de frequências (— 2,0), o diagrama 
é, por razões de simetria, percorrido outra vez 
no mesmo sentido. De acordo com a figura, o 


ponto (1,0) é sempre englobado pelo diagrama 
polar representado. A instabilidade desta monta- 
gem, o que constitui um facto conhecido, encon- 


Fig. X.4 


tra-se assim confirmada através da aplicação do 
critério de Nyquist. 
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A Mobil 
mostra-lhe 


inumeras maneiras 


de gastar 
menos óôleo. 


Reduzir os orçamentos de lubrificação é apenas uma parte 
do Programa Mobil para Redução de Custos. As economias 
totais que se estendem a outros campos podem exprimir-se 
por muitos algarismos. 


Qualquer pessoa pode vender óleos. Mas 
a habilidade estêá em vender menos óleo - e 
nós fizemos disso um sistema: o PROGRA- 
MA MOBIL PARA REDUÇÃO DE CUSTOS. 

Isto significa serviço: serviço técnico, ser- 
viço laboratorial, serviço de treino, serviço 
de informação administrativa... e podíamos 
citar ainda mais vantagens. 

Num cliente importante da indústria por- 
tuguesa conseguimos evitar num só ano o 
consumo de cerca de 200.000300 de óleo 
através de recomendações especiais rela- 
tivas à utilização mais racional dos lubri- 
ficantes. 

Mas poupar óleo foi só o princípio. À 
lubrificação apropriada evitou perdas de 
tempo no uso dos equipamentos e este 


trabalho, mais eficiente, também reduziu 
despesas. 


Na realidade, o que se poupa pode exce- 
der em muito o custo dos lubrificantes. O 
que explica a razão porque vendemos mais 


vendendo menos. 
Mobil 


vendemos mais vendendo menos 


É. D. U. 691.327: 539412: 627.82 
691.327 : 62018 : 627.82 


CARACTERÍSTICAS MECÂNICAS 


E MINEROPETROGRÁFICAS 


DE TAROLOS DE BETÃO 


EXTRAÍDOS DA BARRAGEM DE MIRANDA 


ESTUDO 


DO APARECIMENTO E DA COLMATAGEM 


DUMA FENDA HORIZONTAL 


RESUMO 


Apresentam-se as tensões de rotura à tracção de 205 
provetes e a caracterização mineropetrográjica de 50 lúmi- 
nas delgadas de tarolos de betão extraídos de três furos de 
sondagem efectuados no corpo da barragem e atravessando 
uma fenda horizontal e diversas juntas de betunagem. 

Tentam-se algumas correlações entre as caracteristicas 
obtidas, 

Conclui-se que a fraca resistência de algumas juntas 
de betonagem foi uma das cousas do aparecimento da fenda 
e que a sua colmatagem posterior foi devida à natureza 
encrustante da água da albufeira, comprovada por anúlise 
quémica e pelo estudo das lâminas delgadas, 


1 — INTRODUÇÃO 


O aparecimento duma fenda horizontal, com 
passagem de água, no contraforte 3E da Barragem 
de Miranda, da Hidro-Eléctrica do Douro 
(HIDOURO), levou esta empresa a solicitar ao 
Laboratório Nacional de Engenharia Civil (LNEC) 
a averiguação das causas determinantes da aber- 
tura da referida fenda. 

A HIDOURO procedeu, entretanto, à abertura 
de três furos de sondagem, à cota 510, com vista 
a eliminar possiveis subpressões (fig. 1) e a pos- 
sibilitar a obtenção de tarolos de betão destinados 
a ensaios ulteriores. 

As duas entidades acordaram na realização do 
seguinte programa de trabalhos: 
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Professor do IS 1, Consultor-investigador do LNEC 


SYNOPSIS 


The Authors present the ultimate tensile strength of 205 
samples and the mineralogic-petrographic characterization 
of 50 thin sections from concrete cores obtained in three 
boreholes in the body of the dam through a horizontal crack 
and some construction joints; try some correlations between 
the characteristics obtained; and conclude that poor strength 
at some construction joints was one of the causes of the 
crack, which was aubsequentiy filled up due to the hardness 
of the reservoir water, confirmed by the chemical analysis 
of the water and the observation of the thin sections, 


a) Ensaios de tracção simples de cerca de 200 
tarolos correspondentes a juntas de betona- 
gem e a betão; 


b) Análise mineropetrográfica de 50 lâminas 
delgadas correspondentes a secções de juntas 
de betonagem e de betão; 


c) Análise química da água da albufeira ten- 
dente a determinar a sua agressividade. 


No que se segue, descrevem-se as análises e 
os ensaios realizados e os resultados obtidos. 
Com base nos mesmos tentaram-se algumas cor- 
relações ligando características que pareceram 
mais directamente dependentes e, por fim, formu- 
laram-se, como síntese, algumas conclusões. 


“* Trabalho apresentado em inglês no Simpósio International de Mecânica das Rochas, Madrid, Outubro de 1968, 
e em português nas 2.ºº Jornadas de Engenharia e Arquitectura do Ultramar, Luanda, Junho de 1969. 
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Fig. 1 


2 — CARACTERIZAÇÃO MECÂNICA DO 
BETÃO 


Para se caracterizar mecânicamente o betão e, 
ao mesmo tempo, possibilitar a formulação de 
correlações com os dados obtidos por observação 
microscópica, foram ensaiados, à rotura por 
tracção simples, 205 provetes de modo a obter-se 
uma amostragem tanto quanto possível completa 
ao longo das três sondagens abertas. Atendendo 
a que, para serem comparativos, os valores da 
resistência à tracção se devem referir a prove- 
tes de comprimentos iguais e as amostras dispo- 
níveis apresentarem comprimentos variáveis, clas- 
sificaram-se e cortaram-se os tarolos de sondagem 
de modo a serem definidas as três classes de 
9 cm, 6cm e 3,5 cm. 

Nos topos de cada um dos provetes foram 
coladas peças metálicas apropriadas que ligavam 
cardans especiais (fig. 2). Os provetes, assim pre- 
parados, foram ensaiados numa máquina uni- 
versal de 6 toneladas (fig. 3). 

Os valores das tensões de rotura obtidos, 
relativos a cada sondagem, constam do Quadro 1. 
Neste mesmo quadro estão inscritos também os 
valores das tensões de rotura depois de corrigidos, 
calculados a partir da redução dos valores médios 
das três classes de provetes à classe de valores 
mais baixos que, como seria de esperar, é a de 9 cm. 
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Com o fim de se poderem comparar os valores 
das tensões de rotura obtidos nos ensaios com 
os mesmos valores depois de córrigidos, repre- 
sentaram-se na fig. 4 os respectivos histogramas. 


Fig. 2 


Como no caso presente interessava pôr em 
destaque os valores correspondentes a juntas de 
betonagem, que constituem sempre zonas de 
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Fig. 3 


menor resistência, seleccionaram-se os valores 
mais baixos obtidos nos provetes correspondentes 
a bandas de 0,5 m, centradas nas juntas de 
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Fig. 4 


betonagem, e traçou-se também na fig. 4 o respec- 
tivo histograma para efeitos de comparação. 
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3 — CARACTERIZAÇÃO MINEROPETRO- 
GRÁFICA DO BETÃO 


Foram escolhidos, por exame microscópico, os 
tarolos mais representativos das juntas de beto- 
nagem, a partir das quais se elaboraram 50 
lâminas delgadas que foram estudadas ao micros- 
cópio polarizante. 

Com este estudo pretendeu-se não só caracterizar 
qualitativamente os diferentes troços das sonda- 
gens referidas, como ter conhecimento das quanti- 
dades dos diversos tipos de inerte e argamassa 
(e singularidades) presentes. 

Para a análise qualitativa das lâminas delgadas, 
usou-se a técnica de identificação habitual nos 
estudos de micropetrografia. Efectivamente o 
inerte, no tipo de betão analisado, é fundamental- 
mente constituído por minerais de rochas grani- 
tóides. Quanto à argamassa, houve que proceder 
a determinações de índole mineralógica para a 
caracterizar cabalmente. 

Prestou-se, ainda, especial atenção à presença 
de singularidades relevantes como seja a exis- 
tência de calcite, de microfissuras, de vazios, etc. 

Na análise quantitativa foram seguidos os 
princípios da análise modal tão usada em petro- 
grafia. 

Para se obter o modo, composição real da rocha, 
no caso vertente a composição real do betão, 
usou-se uma platina integradora. Tenha-se em 
atenção que em petrografia ao conceito de modo 
contrapõe-se o de norma, que se refere à compo- 
sição virtual da rocha. Esta composição virtual é 
a que se obtém a partir da análise química da 
rocha, após cálculo petroquímico, mediante o 
qual se determinam os minerais virtuais represen- 
tativos da análise química de que se partiu. 
Estes minerais virtuais podem corresponder a 
minerais reais que existem efectivamente na rocha 
(e contados no modo) ou não existir na rocha, 
mas representar minerais aí presentes. 

Para o cálculo da composição modal do betão 
usou-se uma platina integradora, por meio da 
qual se realizaram séries de trajectos paralelos 
e equidistantes sobre as lâminas delgadas em 
estudo. 

A relação de Delesse, ou de «área-volume» 
em petrogafia estabelece que a razão entre 
a área ocupada por um mineral e a área total da 
amostra é uma estimativa consistente da percen- 
tagem em volume do mineral na rocha. 
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Aspecto geral do betão da sondagem C. D. 1 a05 0,52 m. À branco Mesmo aspecto da fotografia anterior, mas com os nicois cruzados. 
clastos de quartzo, a cinzento escuro e de forma alongada, palheta de 

biotite (inerte). A argamassa isotrópica (pobre ou desprovida de por- 

tlandite). No centro da fotografia, uma cavidade circular preenchida 

por calcite fibrorradiada, forrando interiormente as suas paredes 


Aspecto particular do betão da sondagem C. D.1,aos 7,35 m. A branco, Pormenor da fotografia anterior mostrando a matriz rica de portlandite 
clastos de quartzo, a cinzento, clastos de quartzo e de feldspatos com (anisotrógica) e os vacuolos com fina película interna de (0H)J2 Ca 
geminação polissintética (canto superior esquerdo). À argamassa é rica 

de portlandite. Os vacúolos circulares (a negro) estão vazios e são for- 

rados por portlandite 


Fig. 9 p x 32 Fig. 10 N x 90 
a x 4 enndanom CP MN AS 1 Af , | atonrafia armar Í fntaniraf FLA A tÃ e dao fi ; r várias A 
Aspecto do betão da sondagem CL D, 2 US | 4 mM Na fotografia | menor jd ) gralia anterior. mem id5 vihgs é Ongas vi Jid5 é 
observam-se finas venulas subparaleilas, encaminhando se do canto infe calcite hã plagas de aicite de contomos irregquiare 


rior esquerdo, para o canto superior direito. Tais vénulas são de cal- 
cite e preenchem fissuras que atravessam o merte (quartzo e feldspatos) 


